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ABSTRACT 

The Carnegie Hubble Program (CHP) is designed to improve the extragalactic dis- 
tance scale using data from the post-cryogenic era of Spitzer. The ultimate goal is a 
determination of the Hubble constant to an accuracy of 2%. This paper is the first 
in a series on the Cepheid population of the Large Magellanic Cloud, and focusses on 
the period-luminosity relations (Leavitt laws) that will be used, in conjunction with 
observations of Milky Way Cepheids, to set the slope and zero-point of the Cepheid 
distance scale in the nrid-infrared. To this end, we have obtained uniformly-sampled 
light curves for 85 LMC Cepheids, having periods between 6 and 140 days. Period- 
luminosity and period-color relations are presented in the 3.6 /im and 4.5 pm bands. 

We demonstrate that the 3.6 pm band is a superb distance indicator. The cyclical 
variation of the [3.6]— [4.5] color has been measured for the first time. We attribute the 
amplitude and phase of the color curves to the dissociation and recombination of CO 
molecules in the Cepheid’s atmosphere. The CO affects only the 4.5 pm flux making it 
a potential metallicity indicator. 

Subject headings: Cepheids — distance scale — infrared: stars — Magellanic Clouds 
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Introduction 


In the past decade, t he “factor-of-two cont r oversy” surrounding the value of the Hubble 


con s tant h as bee n resolved ([Freedman et al. 


2001 


Spergel et al.1 120071 : iFreedman fe Madorel l20ld : 


Riess et al. 2011). However, the uncertainty on Ho currently remains at the level of a few percent 
and dominates the error budget when deriving other cosmological parameters such as the equa- 
tion of state parameter, ujq. Reducing the uncertainty on Ho will_enable us to further constrain 
the parameters in the concordance cosmological model — as Hu ( 2005 ) states: “The single most 
important complement to the CMB for measuring the dark energy equation of state at z ^ 0.5 is 
a determination of the Hubble constant to better than a few percent.” 

The Carnegie Hubble Program (CHP) is a “Warm Spitzer” legacy mission (Spitzer Exploration 
Program PID 60010 “The Hubble Constant”, P.I. Freedman). It has the primary goal of reducing 
the systematic uncertainties on the Hubble constant to 3% using Spitzer data alone, and ultimately 


2% with the addition of JWST. The program uses the 3.6 and 4.5 gm bands of IRAC (Fazio et al 


2004) to recalibrate the Cepheid distance scale in the mid-infrared, moving out in to the Hubble 
flow by calibrating and applying a revised mid-IR Tully-Fisher relation. 


A detailed description of the Carnegie Hubble Program is given in Freedman et ah (2011); an 
overview of our strategy follows below. 

The zero-point of the Cepheid period-luminosity (PL) relation (Leavitt law) will be calibrated 
using Galactic Cepheids. This sample contains the 10 Cepheids with trigonometnc_parallaxes mea- 
sured using the Fine Guidance Sensors on the Hubble Space Telescope (Benedict et al. 2007 ). We 
also observed 27 Cepheids within 4 kpc of the Sun — close enough for GAIA parallax measure- 
ments. Seventeen of these are known members of Galactic open clusters for which main sequence 
fitting distances are available. The Milky Way program is described in Monson et al. (2011, in 
preparation). 

The second rung on the distance ladder (the subject of this work) involves the long-period 
(P > 10 days) end of the Leavitt law. We will calibrate the slope and dispersion of this period 
range using a sample from the Large Magellanic Cloud. The LMC is an important stepping stone 
for the extragal actic distance scale; it was used as the base for the Key Project distance scale 
(IFreedman et al.1 |200l|) , but uncertainty in its distance eventually became the dominating system- 
atic. The calibration found here will then be applied to nearby galaxies in the Local Group, forming 
a Cepheid distance ladder based on a single telescope and instrument. Additionally, by including 
galaxies s uc h as M33, the effect of metallicity on the PL slope and zero-point will be studied (e.g., 
Scowcroft et al.l 20091 ) . The calibration will then be applied at a step further out, extending to 
galaxies beyond the Local Group, such as M81 and NGC 4258, the well known maser galaxy. 

Next, we will extend our distance ladder into the Hubble Flow using the Tully-Fisher relation, 
itself re-calibrated in the mid-infrared. The dispersion of this relation is expected to drop at 
IRAC wavelengths, simply due to the decreased sensitivity to inclination-dependent extinction 
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corrections, and will be calibrated using galaxies that contain Cepheids. The distances to over 500 
target galaxies in our program will be measured. Finally, we will complete our distance ladder 
by adding 44 supernova-host field galaxies. These steps will allow us to measure accurate and 
consistent distances from the closest Cepheids to the most distant galaxies. 


This work deals with Cepheids in the Large Magellanic Cloud (LMC). The LM C Cep heid pop- 
ulation has been studied extensively for over half a century. For example, Shapley & Nail (j.1940) 
presented well-sampled light curves for 40 Cepheids in the LMC, andjrtudy of thepopulation con- 
tinues to this day with large-scale projects such as OGLE III ( Soszvnski et al. 2008j). Its proximity 
means we can readily observe Cepheids over a range of periods, from those that are similar to many 
of the Galactic parallax sample (i.e. P < 10 days), to the long period Cepheids with P > 10 days, 
which more generally overlap Cepheids most easily observed in our most distant targets. 


Until relatively recently, the majority of distance measurements have been undertaken at op- 
tical wavelengths (e.g. F01). Cepheid studies at optical wavelengths have their shortcomings, 
the main one being extinction. Although the effect can, to first order, be removed by use of the 
reddening-free Wesenheit index (Madore 1976) this technique still requires prior knowledge of the 
extinction law, as well as an assumption that it is universal. By moving to the mid-infrared, 
reddening and extinction are diminished by around a factor of twenty (jRieke & Lebofsky 1985), 
making their absolu te co ntribution and their uncertainties negligible. In addition to the drop 
in extinction effects, McGonegal et al. (1982) describe two other advantages of the infrared over 
the optical. In the infrared the amplitudes of the Cepheids 1 lightcurves decrease, as does the in- 
trinsic width of the instabili ty str ip, because these wavelengths are less sensitive to temperature 
changes. Indeed McGonegal et ah (119821 ) demonstrated that the width of the near-infrared H band 
period-luminosity relation from random observations is less than the width from time-averaged B— 
ba nd obser vations. Near-infrared observations of LMC Cepheids were more recently obtained by 


Persson et al. (120041 1. By moving to longer wavelengths in the mid-infrared, in combination with 


well phased observations, we can decrease the measured width even further. 


The IRAC imager on Spitzer is a superb instrument for undertaking a recalibration of the 
Cepheid distance scale. Like Hubble , Spitzer has the advantage of operating in a stable environment 
without weather or seeing variations, and with great flexibility in scheduling. As such, we have 
been able to obtain precise and deterministically well-sampled light curves for 85 Cepheids. The 
observations of the LMC Cepheids were one of the first programs undertaken by post-cryogenic 
“Warm Spitzer”. 

The target selection and observations are described in Section [2] and the photometry and 
calibration are discussed in Section [3j Light and color curves are presented for each Cepheid in 
Section [JJ PL relations, including a discussion on their use in determining the tilt of the LMC, are 
given in Section [5j The period-color relation at mean light has been measured for the first time at 
these wavelengths; it is discussed in Section [6l Section [7] provides a summary. 
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Observations 


2.1. Target Selection 


The LMC has a large sample of well-studied C epheids to choose from, as light curves have been 
measured for many thousan ds. The OGLE studies ( Udalski et al.lll992 : Udalski. Kubiak Sz Szvmanski 
1997; Soszynski et ajj 2008), while having a primary goal of observing micro-lensing events, have 
produced extremely large Cepheid catalogs for the L MC and SMC in the optical. For this study, we 
have chosen a subset of the sample observed by Persson et ah (120041 . hereafter P04), who studied 
92 LMC Cepheids in the infrared JHK S bands. These Cepheids also had optical ( UBVRI ) light 
curves available. The majority of them are free from crowding, making them ideal for this study. 
We chose the 85 Cepheids with periods greater than about 6 days. This ensures that the period 
distribution overlaps with both the Galactic sample, as well as the most remote galaxies in our 
program. The Cepheids are distributed over the face of the LMC, allowing any known inclination 
effects (c.f., P04) to be re-measured and accounted for. 


2.2. “Warm Spitzer” Observations 

The data were obtained in programs 61000, 61004, 61005, 61006, and 61007, executed between 
2009 October 3 and 2010 July 18, at the start of the “Warm Spitzer” mission. Each Cepheid was 
observed at 24 epochs. For those Cepheids with P > 12 d, the observations were taken over a 
single cycle, spaced on average by P/24 days. This produced regularly-sampled light curves while 
preserving the high scheduling efficiency of the spacecraft. The shortest-period Cepheids (P < 
12 d) could not be observed in such a tightly-phased manner; instead we scheduled 24 observations 
separated by 12 ± 4 days. The Appendix explains the choice of 24 evenly-spaced epochs. 

Each IRAC observation consisted of ten frames, five in each of the 3.6 /rm and 4.5 ^m bands, 
with a medium-scale, 5-point gaussian dither pattern. (Dithering of exposures is essential for 
accurate photometry.) The frame time was 2.0 s, giving an effective per-pixel integration time of 
1.2 s in each frame. IRAC is known to be susceptible to persistent images (i.e. latent images of 
bright, but often unsaturated, stars can remain in later frames) and this effect is apparently worse 
in Warm Spitzer than in the original mission (see IRAC handbook, version 2 .Cl|, section 7.2.8). 
By using a large scale dither pattern, where subsequent observations of the Cepheid are spaced by 
approximately half the chip, we have mitigated the possibility of our photometry being affected by 
previous observations. 


: http : //irsa. ipac . caltech. edu/data/SPITZER/docs/irac/ iracinstrumenthandbook/ 
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3. Data reduction and Photometry 


These early “Warm Spitzer” data were some of the first to be reduced with the new “Warm 
IRAC” pipeline. The images were processed using the Spitzer Science Center pipeline version 
S18.18.0, which does not include the artifact correction steps present in the cold mission pipelines. 
Therefore, in addition to the normal Spitzer Science Center (SSC) processing, the BCD (Basic Cal- 
ibrated D ata) frames^ were run through the Warm-Mission Column Pulldown Correctol mopex 
(jMakovoz & Kha n 20051 ) was used with the standard overlap and mosaic name lists for the overlap 
correction and mosaicking. 

The most important objective of this work is to produce an accurate calibration of the Cepheid 
period-lumin osity re lation. Our photometry is on the same scale as the standard calibration de- 
scribed in [Reach et al.l ( 200 51). They define the flux scale and zero point such that they correspond 
to a 10-pixel radius aperture, with a sky annulus running from 12 to 20 pixels (in units of native 
IRAC pixels) with pixel-phase and array-location corrections applied as necessary. 

Although the majority of the Cepheids in our sample are isolated, several have close neighbors 
that fall inside the 10-pixel aperture. These stars could not be measured correctly using aper- 
ture photometry. In order to reduce the entire LMC sample in a consistent way, Point Response 
Fu nctio n (PRF) fitting was used in place of aperture photometry, utilizing the SSC APEX software 
([Makovoz & Marleanl 2005). The point response function (PRF) models the detector response to a 
point source, and can be thought of as an over-sampled representation of the point spread function. 
The BCDs were mosaicked to obtain the bad pixel maps, but the PRF fitting was performed on 
each BCD individually using the apex-lframe name list, with two changes. First, rather than using 
the PRF corresponding to the center of the array, the set of 25 spatially-dependent PRFs was 
used for each channel. The IRAC detectors are known to have some spatial response variations; 
using the PRF map ensured that any variation in flux was properly accounted for. Second, the 
normalization radius - the radius to which the fit is calculated - was set to 1000 pixels. As the 
PRF is sampled at 100 times per pixel, the 1000 pixel normalization radius corresponds to the 
calibration aperture radius. For each BCD, all the stars in each image were measured. 


3.1. PRF calibration 

The construction of the IRAC PRF is a complex task requiring a large number of observa- 
tions of a stable star. Warm Spitzer PRF models are not yet available from the Spitzer Science 
Center; however, it is still possible to obtain accurate photometry using the cold PRF, correcting 
as necessary for the predicted differences between the cold and warm models. The most significant 


“http : //ssc . spitzer . caltech . edu/dataanalysistools/tools/contributed/irac/f ixpulldown/ 
3 PRFs are available from http://ssc.spitzer.caltech.edu/irac/calibrationfiles/psfprf/ 
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difference between the two models is the pixel-phase correction (PPG). 

The PPG is necessary given that the response of the IRAC detectors is known to vary on a 
sub-pixel level. This means that the measured flux of an object not only depends on the position 
of the star in the frame, but also on the position of the peak of the PSF within the pixel. The 
fractional position within a pixel is known as the “pixel phase”. In the cryogenic mission the 
pixel-phase effect was around 4% (peak-to-peak), but is approximately double that in the warm 
mission. In addition to the change in magnitude, whereas the cold mission PPC was represented by 
a symmetric Gaussian centered on the peak response position, the warm mission PPC is described 
by a Gaussian, differing in width in the x and y directions. 

As the PRF is an oversampled representation of the pixel’s response to light, it includes the 
cold mission PPC by definition. To correct the warm mission data, the cold PPC must first be 
multiplied out of the APEX photometry results, and the warm correction then divided in. Although 
it is possible that the pixel-phase effect will cancel out when the results from a dither pattern are 
averaged together, in our data the PPC was applied to each BCD before the averaging took place. 


3.2. Photometric stability 

The stability of IRAC is of fundamental importance to our program. We have quantified the 
stability of the instrumental zero-point by measuring constant stars in our Cepheid fields, and 
monitoring their magnitudes over time. The photometry method was identical to that used on 
the Cepheid targets. In Figure CD we show the photometry for constant stars in two fields. Star 
A was chosen in the HV00872 held, at a = 4 h 55 m 17.09 s , <5 = — 67°27 , 53.82” , and Star B was 
in the HV02883 held at a = 4 h 56 m 21.24 s , 5 = — 64°41 , 35.65”. In both cases, the dispersion of 
the measurements is approximately ±0.01 mag, and no systematic trend with time is seen. The 
Cepheids in the helds have very different periods (29.82 d vs. 108.97 d), hence the time baseline 
and the cadence of the observations differ. This confirms both the short and long term stability of 
IRAC. 


3.3. Photometry comparison 


As mentioned above, the IRAC magnitude system is defined by the hux measured within a 
10-pixel aperture. To verify that the photometry for this work is on the calibrated system, the data 
for the Cepheid HV00872 were also processed using the SSC IDL routine "photo -for star. pro (SSC 
2010, private communication). This routine performs aperture photometry on a single star using 
two aperture sizes — a 3 pixel aperture with a 3-7 pixel sky annulus, and the 10 pixel aperture 
with a 12-20 pixel sky annulus used in the calibration definition given by Reach et alj (120051 ). It 


also includes the latest pixel-phase corrections and array-location corrections. 
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The APEX photometry was compared to the results from both aperture sizes. No systematic 
trends with magnitude were found. The difference between the APEX and the 3-pixel aperture 
fluxes were consistent with the published aperture correction for transformations between the 3 
and 10 pixel apertures. Both the APEX and small aperture results were then compared with the 
large aperture values. As was expected, the scatter in the large aperture results was much larger 
than the other two methods, but no systematics were found. Aperture corrections were derived 
for each band, and agreed with the cold mission values (1.124 and 1.127 at 3.6 n m and 4.5 /mi 
respectively) to within 1%. 


To put the PRF photometry on the same system as the standard aperture photometry the 
following correction was used: 


F = 


Fprf x A X P PC Cold 


PPCw 


arm 


(i) 


where F is the flux in the standard system, Fprf is the flux output by apex, A is the aperture 
correction (1.03 and 1.02 for channels 1 and 2 respectively) and PPCcoid and PPCwarm are the 
cold and warm pixel phase corrections. The fluxes were converted to magnitudes using the standard 
zero-magnitude flux densities of 280.9 Jy ([3.6]) and 179.7 Jy ([4.5]). 


Results 


4.1. Periods 


We first refined the available peri ods as follows. Light curves were produced in the 3.6 and 
4.5 gm bands. The data points from iMadore et al.1 (J2009|), taken as part of the SAGE project 
(Meixner et ah 2006), were also phased in. The periods given in P04 were taken as initial values; 
these were improved by iterating on the period and re-phasing the data until a minimum dispersion 
in the light curves was obtained in all available wavebands (U. B,V,R, /, J, H , K. [3.6] , [4.5] , [5.8] , [8.0] , 
but note that not all wavelengths were available for all Cepheids). On average, the periods changed 
by less than 1%. In some cases it was clear that the period of the Cepheid had changed in the 
time interval between the optical and infrared observations. This is to be expected, as we utilized 
as much archival data as possible; the time baseline for some Cepheids was of the order of many 
decades, which is comparable to the time many Cepheids have been found to significantly change 
their periods (for a discussion of period changes in a sample of LMC Cepheids, see Pietrukowicz 
2001). The practice of fitting all wavebands simultaneously makes it easy to find objects for which 
this is the case, as the period that fits the later observations will result in an incorrectly phased 
light curve for the archival data, and vice versa. When this was the case the IRAC and JFLK data 
alone were used to refine the period and to phase the data. The finally adopted periods are given 
in Table 03 


4.2. Infrared Light Curves and Mean Magnitudes 


Once the periods had been adjusted, the IRAC light curves were plotted using just the CHP 
data. This avoids issues such as slight period changes and larger photometric uncertainties in the 
archival data that could affect the IRAC fit. The resulting light curves are shown in Figure [2] in 
order of increasing period. The GLOESS program (described in section 3.2 of P04) was used to 
interpolate between the data points and to thus draw the smooth curves. GLOESS was also used to 
determine tinre-averaged (intensity-mean) magnitudes and uncertainties. The results are given in 
Table |T] 

For the [3.6] (top) and [4.5] (middle) panels, the statistical uncertainties on the individual 
measurements are comparable to the size of the plotted points. Table [2] gives a sample of the 
individual data points for HV00872; all data for this program are available with the online version 
of the article. 


The bottom panels of the plots in Figure [2] show the [3.6] — [4.5] color variation through the 
Cepheid pulsation cycle. As light curves of this quality have not previously been available in 
the near-IR, this cyclical effect has not been noted before. The effect is almost certainly du e to 


absorption by the 4.65 //m CO bandhead in the 4.5 ^m filter, as first suggested by Marengo et al 
(|201 0al). The color variations are discussed below in Section 


5. Mid-infrared Period— Luminosity Relations 

Period-luminosity (PL) relations for the two IRAC_bands were fitted using the weighted-least- 
squares technique described in lAkritas &; Bershadv (119961 ). This takes into account the intrinsic 
width of the PL relation, giving a more robust result than a regular weighted-least-squares fit. 


The PL relations take the form 


M = a(log P — 1) + b 


(2) 


and were fit using four different period cuts. The results of the different fi ts are give n in Table [3j 
The relations derived from the two phase points available from SAGE (M adore et al. 20091 ) are 
also given for comp ariso n. In Table [3j N i s the number of Cepheids in the fit, a and b were both 
found using the Akritas & Bershady (199q) weighted-least-squares method, a a and cp, are their 
respective uncertainties, and a is the standard deviation of the fit. 


The finally adopted relations are plotted in Figure [3] The solid line shows the adopted PL 
relations; these were the fits that used only Cepheids in the range 10 d < P < 60 d. Note the 
systematic deviation of the longest-period Cepheids, which are well known to be significantly fainter 
than the PL relation at virtually all wavelengths. Consistent with other studies (such as F01) we 
have excluded objects with P > 60 days from the fitting. The SAGE PLs (given in Table [3]) are 
entirely consistent with the values found in this work. Although we find no evidence in our data to 
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suggest a break or change in slope in the P = 10 days region (c.f., ISandaee. Tammann fe Reindl 
(2004!))) we have also cut the low end of the period distribution, again consistent with F01. 


In both cases, the dispersion around the relation has fallen dramatically in comparison with 
previously published determinations. The previous measurements from SAGE gave standard devi- 
ations of 0.135 and 0.141 mag in [3.6] and [4.5] respectively, whereas the scatter around_our final 
relations are 0.108 and 0.115 mag. This is because the magnitudes used in the Madore et al. (2009) 
PL relation are not fully phase-averaged, but rather they are the simple mean of two random phase 
points (see the Appendix). 


The precision of our phase-averaged light curves is best demonstrated by comparing the resid- 
uals from the [3.6] PL relation to those from the [4.5] relation. This is shown in Figure |U In this 
figure the small inset plot shows the entire data set, with the main plot showing just the objects 
with A mag < ±0.2, which encompassed 80 of the 85 Cepheids in the sample. It is clear that the 
residuals are highly correlated. That is, if a star lies above the mean PL relation in one band, it 
will tend to be above the mean in the other band as well. The same effect in the J and K$ bands 
was found by P04. 


Excluding two outliers, the range in A mag is ±0.2 mag. This spread is a combination of two 
effects: depth and tilt effects in the structure of the LMC, and the intrinsic magnitude width of 
the instability strip. The residual dispersion in this correlation is only 0.02 mag, suggesting that 
we are seeing the photometric limit of the data. 

Figure [5] shows how the residuals of the ([3.6] — [4.5]) color correlate with the PL residuals 
of the [3.6] magnitude. As we know from the previous plot that any dispersion in Figure [5] arises 
from the intrinsic properties of the Cepheids rather than photometric uncertainties, the width of the 
instability strip can be measured from this figure. Assuming that the instability strip is represented 
by a rectangular distribution, i.e. it is uniformly filled and has hard limits at the blue and red edges, 
the dispersion around the mean is related to the width by the relation a = i?/\/l2, where R is the 
full width of the distribution. The color dispersion of the data is ±0.018 mag, leading to a color 
width of the instability strip of 0.062 mag. 


5.1. Tilt of the LMC 

By combining our knowledge of the Cepheid’s location within the LMC with its position in 
the period-luminosity plane, the tilt of the galaxy with respect to the plane of the sky can be 
constructed, as was done in P04. 

In Figure [ 6 ] we show the positions of the Cepheids within the LMC in projection on the sky. 
The bottom left panel shows the position of each of the Cepheids relative to the center of the LMC. 
Solving for PL coefficients a and b, but now including ^ and 77 (distances from the center of the 
LMC in RA and Dec) as free parameters, we can obtain the tilt and the angle of the line of nodes if 
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the Cepheids on average lie in a plane. The top left and bottom right plots show the distribution of 
the points when viewed from orthogonal directions in the plane of the sky. Empty circles represent 
the Cepheids that are brighter than the average (i.e. they lie above the mean PL relation), and 
filled circles show those that are fainter than the average (below the mean PL relation). In this 
figure the reference coordinate system has been rotated by 50° so that the line of nodes of the LMC 
is parallel with the x-axis. In this plot, it is clear that the Cepheids in the bottom half tend to 
be fainter than the PL, hence further away, while those in the top half tend to be systematically 
brighter, hence closer than average. The simplest interpretation is that the LMC is in fact tilted 
with respect to the plane of the sky. 

By tilting the LMC by 28° relative to the plane of the sky we can minimize the width of the 
histogram in the top right corner of Fig. [6l This is in good agreement with the results of P04, and 
will be the subject of a more complete discussion of the data at a later time. 


6. The Role of CO Absorption 


We believe that the variations in the [3.6] — [4.5] color, both at mean light and throughout 
the pulsation cycle, are due to CO absorption in the 4.5 fim band. We do not see any evidence of 
extended emiss i on sur roundi ng the Cepheid s in ou r images. This is consistent with the results of 
Marengo et al.1 (|2010al lbl) and lBarmbv et al.1 (|201ll ) who found evide nce for circumstellar em ission 
in the longer wavelength IRAC bands, but nothing at 3.6 or 4.5 /mi. Maren go et ah] (2010a) notes 
that the [3.6] — [4.5] color of the Cepheids is indistinguishable from their non-variable control stars, 
and is therefore caused by an intrinsic feature of stars of this spectral type, regardless of whether 
or not they are Cepheids. This rules out circumstellar envelopes and mass-loss as the source of the 
feature, leaving only the CO bandhead. 


6.1. Period— Color Relation 


The possibili ty of a peri od-color (PC) relationship at mid-infrared wavelengths was correctly 
noted by Marengo et al. (2010a). However, their magnitudes (and therefore colors) were measured 
at only a single epoch, introducing large scatter into the PC relation, as is shown in their figure 
5. Figure [7] shows the relationship between period and [3.6] — [4.5] color for our sample. The 
period-color relation was fit using Cepheids with 1 < logP < 1.8, and is found to be 


[3.6] - [4.5] = -0.087(±0.012)(log P - 1) + 0.005(±0.005). 


( 3 ) 


The LMC Cepheids clearly exhibit a period-color relation in this period range. The standard 
deviation around the relation is ±0.018 mag. The relation appears to invert for the longest — 
period (log P > 2) Cepheids. Ignoring the log P > 2 sample, the long period Cepheids are brighter 
and bluer in [3.6] — [4.5] than their short period counterparts. At longer periods the Cepheids are 
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at cooler average temperatures, and so the carbon and oxygen are situated in CO molecules. As 
discussed previously, the CO absorbs some of the 4.5 /rrn flux, making the star appear bluer. As 
we move to shorter periods the average temperature rises, the CO is gradually dissociated, and 
the average colors become redder. The period-color relation reverses at P > 100 days. We do not 
know what causes this reversal, but the Cepheids at these periods are deviant in their luminosities 
and colors at almost all wavelengths. 


6.2. [3.6] — [4.5] Color Curves 

The period-color relation is closely related to the variations in color shown in Figure El as they 
are both due to CO affecting the 4.5 /irn flux. In a static stellar atmosphere, with a temperature 
and surface gravity corresponding to an F-type supergiant, the 4.5 /im flux is suppressed. In the 
case of Cepheids, however, we are not dealing with static atmospheres. As the Cepheid contracts 
its temperature increases. As the star becomes hotter, CO in the atmosphere dissociates, causing 
less absorption at 4.5 /im, making the star then appear redder than average. Once the Cepheid 
begins to expand again its temperature drops, allowing the C and O to recombine. As more CO 
is present, the 4.5 jim flux is once again suppressed, making the star bluer, in the particular color 
combination ([3.6] — [4.5]) discussed here. Inspection of Figure [2] shows that the phasing of the light 
and color curves exhibits precisely this behavior. 

There also appears to be a relation between period and the amplitude of the color variation, 
such that the short period Cepheids (P < 10 days) tend to have flat color curves. The cyclical 
variation appears to “turn on” around P = 10 days, increasing in amplitude with period. Again, 
this is a temperature effect. The short period Cepheids are so hot that the CO is always dissociated, 
whereas the longest period Cepheids have lower maximum temperatures and so never dissociate all 
their CO, thus have the largest color amplitude. 

The dependence of color on CO implies that the 4.5 fi m band should be subject to a metallicity 
effect, effectively ruling it out as a distance indicator. Fortunately, the bandhead is confined to 
the [4.5] filter alone; no CO or other molecular feature is present at 3.6 /im, and so it can still be 
confidently used to measure Cepheid distances. The physics behind the cyclical CO variation, along 
with the effect of metallicity on Cepheid’s IRAC colors will be described in detail in Scowcroft et 
al. (2012, in preparation). 


6.3. PL Slopes and CO Absorption 


We now return to consideration of the slopes of the PL at different wavelengths. Figure [8] shows 


the dependence of the slope on wavelength, a plot that is repeated from figure 4 of Freedman et al 


12003). The [3.6] and [4.5] slope values have been updated using the Table 3 values in this paper. 
The slope of the [4.5] PL relation is clearly deviant. We can again understand this as due to CO 
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as follows: as we move from shorter to longer periods we are moving to cooler stars; these have 
more CO absorption on average. More CO means depressed 4.5 /xm flux and thus a shallower PL 
slope. Therefore, this plot gives further evidence that the color curves are due to CO absorption 
that affects only [4.5]; the [3.6] magnitude is unaffected. 

A smaller magnitude effect is also seen in the K band where a shallower CO bandhead is 
present. However, this effect can be mitigated using the K$ band, which excludes the longest 
wavelength section of the bandpass — the part that happens to include the CO bandhead. 


7. Summary 

We have used the IRAC camera on “Warm Spitzer” to measure [3.6] and [4.5] light curves for 
85 Cepheids spread across the face of the LMC and covering a 6 to 140 day period range. Each 
star has been measured at 24 phase points. Our data agree well with the random phase data of the 
SAGE project, the latter taken for different purposes. 

The light curves (shown in Figure [2]) allow us to measure the mean magnitudes and color of each 
Cepheid to a precision never achieved before in the mid-infrared. These lead to our best PL relations 
based on the measurements of the [3.6] and [4.5] magnitudes for 67 stars with 10 < P < 60 days. 
The PL relations have been measured to such high precision that their scatter directly correlates 
with the front-to-back geometry of the sample of Cepheids in the LMC. This correlation recalls 
a similar effect presented in P04. Using the residuals of the period-luminosity relation, we have 
measured the tilt of the galaxy as 28°, with respect to the plane of the sky. 

The dense and uniform sampling of our Cepheid light curves reveals a cyclical color variation 
which we interpret as being due to the effect of temperature-sensitive molecular absorption of CO 
in the Cepheid’s atmosphere. When the Cepheid is cool, the CO present in its atmosphere absorbs 
at 4.67 jum, suppressing flux exclusively in the IRAC 4.5 /jrn band, making the star appear blue 
in the [3.6] — [4.5] color. As the Cepheid contracts, and so heats up, the CO molecules begin to 
dissociate, ending the suppression of the 4.5 /zm flux and reddening the color. The temperature 
dependence of the CO absorption also explains the period-color relation we observe (see Figure E]) : 
the cooler, longer period Cepheids appear bluer (more highly absorbed by CO at 4.5 /.tin) than their 
short period counterparts. The implications of this temperature-sensitive effect will be studied in 
detail in Scowcroft et al. (2012, in preparation). 

The period-luminosity relations measured here will be used in combination with those from 
our Galactic sample (Monson et al. 2011, in preparation), to establish the bottom rung of the 
extragalactic distance ladder. They will then be combined with parallel observations of Cepheids in 
the Local Group and beyond, to establish the distance scale to a systematic accuracy not achievable 
before now. We will also be testing for the effect of metallicity on the mid-IR period-color relation, 
by combining this sample with Galactic and SMC Cepheids, (Scowcroft et al. 2012, in preparation). 
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9. Appendix 


The observing plan we adopted using Spitzer to optimally follow up known Cepheids is very 
differently constrained from what one would or could attempt from the ground. Ground-based 
observations are obviously restricted by 24-hour resonances with the day-night observing cycle. 
They are often additionally modulated by the lunar cycle (especially for optical observations). 
They are stochastically interrupted by weather, instrument availability, technical problems, and 
time assignments. Randomly spaced observations of known Cepheids may thus be the only way to 
fully sample light curves and to determine their mean values. In the case of space-based observations 
of known Cep heids randomly samp led data may be the only kind av ailable (e.g. data from the 
SAGE project (IMeixner et al.ll2006l h which was mined for Cepheids bv lFreedman et al.1 (|20081 ) and 
Madore et al. (2009)). 


However, the 24-hour availability of a solar-orbiting telescope opens up a unique opportunity to 
design the observing program to be optimally scheduled on a star-by-star basis. Let us first consider 
the case of random sampling. In randomly sampling any given distribution function the error on the 
derived mean value goes down in inverse proportion to the square root of the number N of randomly 
timed observations (i.e., a oc 1/y/lV). What is not widely appreciated is that convergence of the 
error on the mean can be very significantly accelerated by deterministic (non-random) sampling. In 
the case of an individual Cepheid where the time distribution of luminosities (i.e., its amplitude) is 
bounded (but not a priori known), and where th e cycl ing time (i.e. . it s pe riod, but not necessarily 
its phase) is also known, then it is possible (see iMado re & Freedman 2005 ) to construct a simple 
observing strategy where the errors on the mean decrease in direct inverse proportion with the 
number of observations (i.e. a oc 1 /N). Not surprisingly, for a specified number of observations N, 
and a known period of variation P, the optimal separation of data points in time is the uniform 
distribution P/(N+1). 


The distribution function in this case is closely approximated by a rectangular probability 
density distribution. The equivalent sigma for such a distribution is A/\/l2, where A is the total 
range of the distribution function; in the case of a Cepheid, A is its amplitude. In the rnid-IR the 
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amplitudes of Cepheids can reach 0.8 mag (e.g. HV00883), which converts to an equivalent sigma 
of ±0.23 mag. If the required precision on a determination of the mean magnitude were to be set 
at 0.01 mag (i.e. , one percent) then random sampling of the light curve would require over 500 
observations! To reach the same precision using deterministic (uniform) sampling a grand total of 
23 observations are required. To compensate for slight variations away from exact scheduling we 
requested and obtained 24 observations per LMC Cepheid. The light curves shown in Figure [2] are 
a testament to the care and efficiency with which Spitzer was scheduled for this program. 
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0 20 40 60 80 100 


MJD - 55133 


Fig. 1. — Photometry of two randomly chosen constant stars in the HV00872 (Star A) and HV02883 
(Star B) fields. The dashed lines show the average magnitude of the star and ±1 a limits. The 
dispersion in both cases is approximately 0.01 mag. The time baselines differ as the Cepheids in 
the fields were observed with different cadences. No systematic trends are seen. 
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HV00914 P = 6.879 



HV13048 P = 6.853 



HV02337 P = 6.863 




Fig. 2. — IRAC light curves of the LMC Cepheid sample, in order of decreasing period. Point sizes 
are comparable to the uncertainties in the [3.6] (top) and [4.5] (middle) panels. The bottom panels 
show the variation of the IRAC color with phase. Upwards (more negative) in color corresponds 
to greater CO absorption as discussed in Section [HI 
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log P (days) 


Fig. 3. — LMC mid-infrared PL relations. The solid lines indicate the relations given in Table 02 
cut at 10 and 60 days, i.e. the open circles are not included in the fit. The dashed lines show ±2 a. 
Error bars are consistent with the size of the points. 
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Fig. 4. — Magnitude residuals from the [3.6] and [4.5] period-luminosity relations, with short period 
(P < 10 days) Cepheids as red open stars, long period (P > 100 days) as blue open triangles, and 
the PL fitting sample as black filled circles. The inset plot shows the whole sample; the main plot is 
cut at ±0.2 mag to show an expanded view of the details. The dashed line indicates A [3.6] = A [4.5], 
the solid line is a fit to the black points. The deviation from a one-to-one fit is believed to be due 
to CO pushing the [4.5] points further from the adopted PL. 
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A([3.6] - [4.5]) 


Fig. 5. — Color and magnitude residuals from the period luminosity relation. The dashed lines show 
the lcr dispersion. Excluding the longest period (logP > 2) stars, the color dispersion around the 
[3.6] ridge line is 0.018 mag. Assuming that the instability strip has a regularly filled, rectangular 
distribution, the width will be vT2o\ This gives the color width of the instability strip as 0.062 


mag. 
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TILT OF THE LMC 



* 

Fig. 6. — Positions of the Cepheids within the LMC. Open circles show Cepheids which lie above 
the PL relation (i.e. are brighter than average), filled circles show Cepheids which lie below it. 
The bottom left plot shows their positions relative to the center of the LMC, with the coordinate 
system rotated such that the line of nodes is horizontal. The top left and bottom right plots show 
the deviations from the PL relation separated into the £ and ij components. The top right plot 
shows the histogram of the size of the PL deviations. The dispersion in the PL is minimized by 
assuming a tilt of 28° in r/ relative to the plane of the sky. 
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Fig. 7. — Period-color relation for our LMC sample. The solid line shows the PC relation of 
Equation^ with the dashed lines showing Tier. The trend apparently reverses for the very longest 
period (log P > 2) objects. 
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Fig. 8. — Relationship between effecti ve wav elength and the slope of the Cepheid PL relation. 
The PL sl opes are taken from Madore &; Freed m a, n (1998) ( BVRI ), P04 ( JHK ), this work ([3.6], 
[4.5]) and (Mador e et al, (20 09) ([5.8], [8.0]). There is an asymptotic trend towards —3.45, which is 
consistent with the theoretical prediction of a slope of —3.4 for the period-radius relation. However, 
the [4.5] (and to a smaller extent, K ) slope is clearly discrepant. We attribute this to the CO 
bandhead that suppresses flux in the 4.5 n m. The [3.6] slope is not affected — it follows the trend 
of the other wavelengths — and is therefore appropriate to use for distance determinations. 
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Table 1. Magnitudes for LMC Cepheids. 
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Table 1 — Continued 


Cepheid 

Period 

(days) 

[3.6] 

(mag) 

0"[3.6] a 

[4.5] 

(mag) 

014.5] a 

[3.6] - [4.5] 
(mag) 

cr [3.6]— [4.5] a 

HV01005 

18.716 

11.922 

0.007 

11.897 

0.006 

0.023 


HV01006 

14.217 

12.162 

0.005 

12.172 

0.005 

-0.010 


HV01013 

24.131 

11.392 

0.006 

11.458 

0.005 

-0.064 


HV01019 

13.659 

12.252 

0.003 

12.253 

0.003 

0.001 


HV01023 

26.559 

11.386 

0.006 

11.441 

0.006 

-0.052 


HV02244 

13.977 

12.219 

0.004 

12.230 

0.005 

-0.011 


HV02251 

27.887 

11.124 

0.007 

11.166 

0.007 

-0.043 


HV02257 

39.395 

10.722 

0.006 

10.784 

0.006 

-0.062 


HV02260 

12.987 

12.633 

0.005 

12.620 

0.005 

0.012 


HV02262 

15.832 

12.035 

0.004 

12.067 

0.004 

-0.034 


HV02270 

13.625 

12.399 

0.004 

12.424 

0.004 

-0.025 


HV02279 

6.895 

13.077 

0.013 

13.050 

0.014 

0.023 


HV02282 

14.677 

12.218 

0.005 

12.221 

0.005 

-0.003 


HV02291 

22.317 

11.667 

0.006 

11.706 

0.006 

-0.036 


HV02294 

36.552 

10.663 

0.006 

10.700 

0.006 

-0.036 


HV02324 

14.466 

12.205 

0.005 

12.221 

0.004 

-0.017 


HV02337 

6.863 

13.247 

0.016 

13.212 

0.015 

0.034 


HV02338 

42.200 

10.580 

0.007 

10.640 

0.007 

-0.057 


HV02339 

13.880 

12.139 

0.004 

12.143 

0.004 

-0.005 


HV02352 

13.632 

12.282 

0.003 

12.267 

0.003 

0.014 


HV02369 

48.377 

10.297 

0.006 

10.355 

0.006 

-0.058 


HV02405 

6.924 

13.309 

0.011 

13.283 

0.012 

0.025 


HV02432 

10.918 

12.446 

0.015 

12.427 

0.016 

0.018 


HV02447 

118.350 

9.440 

0.006 

9.438 

0.005 

0.003 


HV02463 

13.965 

12.142 

0.004 

12.134 

0.004 

0.006 


HV02527 

12.950 

12.452 

0.005 

12.456 

0.005 

-0.004 


HV02538 

13.869 

12.246 

0.003 

12.266 

0.003 

-0.020 


HV02549 

16.220 

11.885 

0.004 

11.891 

0.004 

-0.006 


HV02579 

13.428 

12.121 

0.003 

12.114 

0.004 

0.006 


HV02580 

16.918 

11.879 

0.006 

11.910 

0.005 

-0.031 


HV02733 

8.722 

12.887 

0.011 

12.855 

0.010 

0.029 


HV02749 

23.105 

11.625 

0.004 

11.682 

0.004 

-0.058 
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Table 1 — Continued 


Cepheid Period [3.6] cr [3 . 6] a [4.5] a [4 . 5 f [3.6] - [4.5] o- [3 . 6 ]-[4.5] a 

(days) (mag) (mag) (mag) 


HV02793 

19.220 

11.734 

0.006 

11.767 

0.005 

-0.033 


2 

HV02827 

78.802 

9.744 

0.004 

9.824 

0.003 

-0.081 


2 

HV02836 

17.528 

11.988 

0.005 

12.008 

0.005 

-0.019 


2 

HV02854 

8.635 

12.810 

0.016 

12.785 

0.013 

0.024 


5 

HV02883 

108.970 

9.953 

0.010 

9.958 

0.009 

-0.006 


2 

HV05497 

99.100 

9.345 

0.004 

9.376 

0.003 

-0.030 


2 

HV05655 

14.211 

12.308 

0.005 

12.325 

0.005 

-0.015 


2 

HV06065 

6.838 

13.301 

0.011 

13.282 

0.013 

0.014 


5 

HV06098 

24.237 

11.224 

0.003 

11.218 

0.003 

0.006 


1 

HV08036 

28.369 

11.329 

0.006 

11.374 

0.006 

-0.044 


2 

HV12452 

8.736 

12.837 

0.016 

12.832 

0.015 

0.003 


5 

HV12471 

15.863 

12.178 

0.004 

12.208 

0.004 

-0.028 


2 

HV12505 

14.389 

12.327 

0.005 

12.338 

0.005 

-0.013 


1 

HV12656 

13.400 

12.269 

0.002 

12.248 

0.002 

0.020 


1 

HV12700 

8.153 

12.975 

0.011 

12.953 

0.012 

0.022 


5 

HV12717 

8.843 

12.877 

0.015 

12.842 

0.015 

0.034 


5 

HV12724 

13.744 

12.419 

0.004 

12.444 

0.005 

-0.025 


2 

HV12815 

26.135 

11.212 

0.006 

11.257 

0.006 

-0.042 


2 

HV12816 

9.108 

12.890 

0.011 

12.869 

0.011 

0.019 


2 

HV13048 

6.853 

13.165 

0.011 

13.152 

0.013 

0.010 


5 

P71-U01 

22.558 

11.661 

0.007 

11.702 

0.007 

-0.038 


2 


a For Cepheids with P < 12 days the uncertainties scale with 1 / \fN rather than 
1/N as the observations of short period Cepheids were not phase locked. 
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Table 2. Photometry of HV00872. 


MJD b 

[3.6] b 

(mag) 

°l3.6] 

[4.5] b 

(mag) 

o'p.s] 

55133.665 

11.538 

0.013 

11.633 

0.015 

55134.756 

11.568 

0.013 

11.605 

0.014 

55136.253 

11.510 

0.013 

11.498 

0.013 

55136.953 

11.485 

0.012 

11.471 

0.014 

55138.621 

11.396 

0.013 

11.372 

0.013 

55139.561 

11.378 

0.012 

11.330 

0.013 

55140.963 

11.334 

0.011 

11.296 

0.013 

55141.986 

11.268 

0.011 

11.276 

0.012 

55143.049 

11.244 

0.012 

11.255 

0.012 

55144.443 

11.205 

0.011 

11.235 

0.013 

55146.159 

11.182 

0.011 

11.232 

0.012 

55146.916 

11.162 

0.011 

11.235 

0.012 

55148.340 

11.165 

0.011 

11.243 

0.012 

55149.495 

11.167 

0.011 

11.267 

0.014 

55151.036 

11.171 

0.011 

11.290 

0.012 

55151.989 

11.187 

0.011 

11.292 

0.012 

55153.002 

11.220 

0.011 

11.315 

0.013 

55155.321 

11.266 

0.011 

11.394 

0.013 

55155.906 

11.282 

0.011 

11.378 

0.013 

55156.963 

11.302 

0.011 

11.419 

0.013 

55158.231 

11.357 

0.012 

11.460 

0.013 

55159.318 

11.389 

0.012 

11.528 

0.014 

55160.793 

11.439 

0.012 

11.523 

0.014 

55161.870 

11.472 

0.012 

11.584 

0.014 


Note. - The Cepheid photometry is avail- 
able in the electronic edition of the Astrophys- 
ical Journal. 

a MJD = JD - 2,400,000.5 d 

b Each magnitude point corresponds to a 
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flux weighted average of the five dither po- 
sitions. 



Table 3: Period-luminosity relations 


Band 

Sample 

N 

a 

CTa 

b 


a 

[3.6] 

lOd < P < 60d 

67 

-3.41 

0.08 

12.74 

0.03 

0.108 


SAGE a 

70 

-3.40 

0.07 

12.73 

0.02 

0.135 


No cut 

85 

-3.16 

0.05 

12.67 

0.02 

0.134 


6 d < P < 60 d 

82 

-3.31 

0.05 

12.74 

0.01 

0.105 


P > 10 d 

72 

-3.14 

0.07 

12.67 

0.03 

0.142 

[4.5] 

lOd < P < 60d 

67 

-3.31 

0.08 

12.74 

0.03 

0.115 


SAGE a 

70 

-3.35 

0.07 

12.77 

0.03 

0.141 


No cut 

85 

-3.10 

0.05 

12.67 

0.02 

0.131 


6 d < P < 60 d 

82 

-3.22 

0.05 

12.70 

0.02 

0.112 


P > 10 d 

72 

-3.11 

0.07 

12.68 

0.03 

0.139 


The SAGE data from Madore et al. ( 20091 ) were re-fit to be consistent with the CHP analysis. 


